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The  objective  of  this  investigation  was  to  continue  the  assess¬ 
ment  of  the  Calspan  perforated-wall,  segmented-plenum  adaptive-wall 
wind  tunnel  for  flows  in  which  the  tunnel  walls  are  supercritical. 

An  important  aspect  of  the  investigation  was  the  development  of  a 
static-pipe  measuring  technique  for  determining  the  normal  velocity 
distributions.  This  technique  was  developed  to  overcome  the  limita¬ 
tions  on  the  number  of  flow-angle  probe  measurements  that  could  be 
made  in  the  Calspan  tunnel.  The  flow  about  a  static  pipe  in  the 
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20.  ABSTRACT,  Concluded. 

presence  of  model/wall-lnduced  disturbances  was  analyzed  within  the 
framework  of  slender-body  theory  to  give  the  relationship  between 
differential  pressure  measurements  across  the  pipe  and  the  steam- 
wise  derivative  of  the  normal  velocity  at  the  pipe  centerline. 
Experiments  with  a  4%  blockage  NACA  0012  airfoil  model  at  Ivbo  •=  0.8 
and  a  "  4°  demonstrated  that  suitably  constructed  static  pipes,  in 
conjunction  with  several  flow-angle  probes,  provide  an  improved 
technique  for  measuring  the  normal  velocity  distributions.  Two 
converging  iterative  steps  toward  unconfined  flow  were  accomplished 
at  this  test  condition  and  established  a  basis  for  plenum-pressure 
control  valve  settings  at  thigher  Mach  numbers.  In  order  to  exam¬ 
ine  a  case  in  which  shock  waves  extend  to  the  walls,  the  Mach  num¬ 
ber  was  increased  to  “•  0.9,  with  the  model  still  at  a  =  4°.  In 
the  initial  tests  at  this  condition,  the  flow  in  the  test  section 
was  choked  downstream  of  the  model.  However,  with  Increased  suc¬ 
tion  Just  upstream  of  the  model,  the  choking  was  relieved  and  ad¬ 
justment  of  the  flow  to  a  reasonable  first  interative  step  was  ac¬ 
complished.  Sufficient  control  appears  to  be  available  to  continue 
interation  toward  unconfined  flow  for  this  case. 
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PREFACE 


The  work  reported  herein  wae  conducted  by  Calspan  Corporation  for  the 
Department  of  the  Air  Force  at  the  request  of  Headquarters,  Arnold  Engineering 
Development  Center /DOT,  Arnold  Air  Force  Station,  Tennessee,  under  Contract 
Number  F40600-78-C-0003 .  Air  Force  technical  representatives  were  Mr.  E.  R. 
Thompson  and  Mr.  A.  E.  Money.  The  Investigation  was  conducted  In  the  Aerodynamic 

Research  Department  of  the  Calspan  Advanced  Technology  Center,  Buffalo,  New  York, 

under  the  management  of  Dr.  John  C.  Erickson,  Jr. 

The  reproduclbles  used  for  this  report  were  supplied  by  the  authors, 

J.  C.  Erickson,  Jr.,  C.  E.  Hlttllff  and  C.  C.  Daughtry.  Appendix  A,1  was 
prepared  by  Dr.  J.  P.  Nennl  of  Calspan  ATC  on  the  basis  of  his  static-pipe 

analysis  by  matched  asymptotic  expansions.  The  authors  wish  to  acknowledge 

the  contributions  by  Professor  W.  R.  Sears  of  the  University  of  Arizona,  a 
consultant  to  Calspan,  as  well  as  by  Calspan  personnel  Dr.  A.  BU.tter,  J.  Nemeth, 
R.  £.  Phlbbs  and  the  late  A.  F.  Gretch. 
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1.0  INTRODUCTION 


The  influence  of  wall -interference  effects  on  models  tested  in  tran¬ 
sonic  wind  tunnels  has  been  an  open  question  for  many  years,  and  it  is  generally 
acknowledged  that  many  questions  remain  unanswered  today.  The  accepted  pro¬ 
cedure  has  been  to  restrict  the  model  size  to  be  tested  and  to  neglect  the 
presumably  small  residual  interference.  However,  recent  requirements  have 
led  to  testing  larger  models  at  high  lift,  while  stressing  extreme  data  ac¬ 
curacy.  This  gave  renewed  interest  to  methods  for  treating  the  wall- inter¬ 
ference  problem. 

This  state  of  affairs  has  led  to  the  concept  of  an  adaptive-wall 
wind  tunnel  in  which  the  flow  in  the  vicinity  of  the  test-section  walls  is 
actively  controlled  in  order  to  minimize  or  eliminate  wall  interference.  The 
distributions  of  the  disturbance  velocity  components  are  measured  at  discrete 
points  along  an  imaginary  control  surface  in  the  flow  field  near  the  walls, 
but  away  from  the  model.  A  theoretical  formulation  for  the  flow  exterior  to 
this  control  surface,  including  the  unconfined- flow  boundary  condition  that 
all  disturbances  vanish  at  infinity,  is  used  to  establish  the  functional  re¬ 
lationships  which  must  be  satisfied  at  the  control  surface  by  the  measured 
disturbance  velocities.  If  the  measured  velocities  do  not  satisfy  these  re¬ 
lationships,  an  iterative  procedure  provides  a  new  approximation  for  the  -flow 
field  at  the  surface,  and  the  flow  through  the  tunnel  walls  is  readjusted  until 
the  measured  qiiantities  satisfy  the  functional  relationships  for  unconfined 
flow.  In  this  way,  the  best  features  of  theory  and  experiment  are  combined  to 
eliminate  wall  interference. 

A  program  of  research  is  in  progress  at  Calspan  Corporation  to  de¬ 
velop  and  demonstrate  the  feasibility  of  the  adaptive-wall  concept.  This 

1  2 
concept  was  described  in  general  and  was  demonstrated  theoretically  in 

1.  Sears,  W. R.  "Self-Correcting  Wind  Tunnels"  (The  Sixteenth  Lanchester 
Memorial  Lecture.)  The  Aeronautical  Journal,  Vol.  78,  No.  758/759, 
February/March  1974,  pp.  80-89. 

2.  Erickson,  J.C. ,  Jr.,  and  Nenni,  J.P.  "A  Numerical  Demonstration  of  the 
Establishment  of  Unconfined- Flow  Conditions  in  a  Self-Correcting  Wind 
Tunnel"  Calspan  Report  No.  RK-5070-A-1,  November  1973. 
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two-dimensional  incompressible  flow.  A  two-dimensional  adaptive-wall  test 
section  was  designed,  fabricated  and  installed  in  the  circuit  of  the  Calspan 
One-Foot  Transonic  Wind  Tunnel.  The  test  section'^*  consists  of  perforated 
upper  and  lower  walls  with  segmented  plenum  chambers,  each  of  which  has  an 
individual  pressure  control.  A  model  with  an  NACA  0012  airfoil  section  and  a 

C 

6-inch  chord  was  fabricated  and  tested  in  the  Calspan  Eight-Foot  Transonic 

Wind  Tunnel  to  establish  the  airfoil  characteristics  in  unconfined  flow  at  a 

Reynolds  number  of  1.0  x  10^  based  on  chord  length.  Data  were  obtained  over  a 

Mach  number  range  from  0.40  to  0.95  for  an  angle-o£-attack  range  from  -2®  to  8®. 

An  experimental  demonstration  of  the  tunnel  in  adaptive-wall  operation  was 
3-4 

accomplished  for  flows  which  were  supercritical  at  the  model  but  subcritical 
at  the  control  surfaces  and  walls. 


Upon  completion  of  those  experiments,  research  was  initiated  at  Calspan, 

with  AEDC  sponsorship,^  to  extend  the  experiments  to  conditions  where  there  is 

supercritical  flow  at  the  control  surfaces  and  at  the  tunnel  walls.  The  present 

investigation  continues  the  research  along  these  lines  and  comprises  one  aspect 

of  an  overall  theoretical  and  experimental  AEDC  effort  on  transonic  adaptive- 

wall  wind  tunnel  development.  The  overall  effort  includes  theoretical  studies 

7  8 

of  the  convergence  aspects  of  adaptive-wall  wind  tunnels,  ’  two-dimensional 

3.  Vidal,  R. J. ,  Erickson,  J.C.,  Jr.,  and  Gatlin,  P.A.  "Experiments  With 
a  Self-Correcting  Wind  Tunnel"  AGARD  CP  No.  174  on  Wind  Tunnel  Design 
and  Testing  Techniques,  October  1975. 

4.  Sears,  W.R,,  Vidal,  R.J.,  Erickson,  J.C.,  Jr.,  and  Ritter,  A.  "Inter¬ 
ference-Free  Wind-Tunnel  Flows  by  Adaptive-Wall  Technology"  Journal  of 
Aircraft,  Vol.  14,  No.  11,  November  1977,  pp.  1042-1050. 

5.  Vidal,  R.J. ,  Catlin,  P.A.  and  Chudyk,  D.W.  "Two-Dimensional  Subsonic 
Experiments  with  an  NACA  0012  Airfoil"  Calspan  Report  No.  RK-5070-A-3, 
December  1973. 

6.  Vidal,  R.J.,  and  Erickson,  J.C.,  Jr.  "Research  on  Adaptive  Wall  Wind 
Tunnels",  AEDC  Report  No.  AEDC-TR-78-36,  November  1978. 

7.  Lo,  C.F.  and  Kraft,  E.M.  "Convergence  of  the  Adaptive -Wall  Wind  Tunnel" 

AIAA  Journal,  Vol.  16,  No.  1,  January  1978,  pp.  67-72. 

8.  Lo,  C.F.  and  Sickles,  W.L.  "Analytic  and  Numerical  Investigation  of  the 
Convergence  of  the  Adaptive- Wall  Concept"  AEDC  Report  No.  AEDC -TR- 79-55, 
November  1979. 
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adaptive- wall  experiments  in  the  Aerodynamic  Wind  Tunnel  (IT)  at  AEDC  to  reduce 

g 

wall  interference  by  means  of  tvro  different  porous  wall  configurations,  and  an 
application  of  the  adaptive-wall  concept  in  three-dimensional  experiments  with 
a  swept- wing/body/horizontal -tail  configuration  in  the  Aerodynamic  Wind  Tunnel 
(4T)  at  AEDC  wherein  wall  porosity  is  constant  over  each  wall,  but  can  vary  from 
wall  to  wall.^^  An  additional  aspect  of  the  overall  AEDC  program  is  the  develop¬ 
ment  of  a  computer  code  to  model  the  inviscid,  full-potential  flow  over  three- 
dimensional  swept  wings  in  an  adaptive-wall  transonic  wind  tunnel. 

The  results  of  the  initial  Cal span  investigation  for  AEDC  are  described 
in  detail  in  Reference  6.  However,  to  provide  continuity  with  the  present  study, 
the  results  of  that  program  are  summarized  briefly  here.  The  initial  experiments 
with  the  6-inch  (6%-blockage3  model  were  devoted  to  determining  a  practical  mode 
of  operation  when  shock  waves  from  the  model  extend  to  the  walls.  This  was 
necessary  because  the  previously  established  operational  mode  for  subcritical 
walls  produced  choked  flow  at  and  downstream  of  the  model.  The  most  practical 
mode  of  operation  which  resulted  from  the  initial  experiments  consisted  of  two 
steps.  First,  wall  control  was  used  to  obtain  the  desired  distributions  of  the 
streamwise  velocity  component  for  a  subcritical-wall  case.  Then,  the  Mach  number 
was  increased  and  the  wall  control  was  readjusted,  sequentially,  until  the  desired 
supercritical -wall  test  condition  was  achieved.  At  the  high  Mach  numbers  of 
interest,  however,  the  available  wall  control  was  limited  locally,  and  changes 
were  required  in  the  experimental  apparatus. 


9.  Kraft,  E.M.  and  Parker,  R.L,,  Jr.  "Experiments  for  the  Reduction  of  Wind 
Tunnel  Wall  Interference  by  Adaptive-Wal 1  Technology"  AEDC  Report  No. 
AEDC-TR-79-S1 ,  October  1979. 

10.  Parker,  R.L.,  Jr,  and  Sickles,  W.L.  "Application  of  Adaptive  Wall  Tech¬ 
niques  in  a  Three-Dimensional  Wind  Tunnel  with  Variable  Wall  Porosity" 
AIAA  Paper  No.  80-0157,  January  14-16,  1980. 

11,  Mercer,  J,E.,  Geller,  E.W.,  Johnson,  M.L.  and  Jameson,  A,  "A  Computer 
Code  to  Model  Swept  Wings  in  an  Adaptive  Wall  Transonic  Wind  Tunnel" 

AIAA  Paper  No.  80-0156,  January  14-16,  1980. 
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A  method  was  developed  for  analyzing  and  estimating  flow-control 
requirements  for  adaptive-wall  wind  tunnels  with  perforated  walls  and  seg¬ 
mented  plenum  chambers.  The  analysis  treats  the  auxiliary  conqjressor  circuit 
for  wall  control  and  provides  an  approximate  method  for  examining  the  trade¬ 
offs  between  the  compression  ratio  of  the  auxiliary  compressor,  wall  open-area 
ratio,  and  model  size.  For  the  purposes  of  the  previous  study,  the  most 
expeditious  alternative  for  extending  the  operating  range  of  the  Calspan  tun¬ 
nel  was  to  reduce  model  size  by  fabricating  a  4-inch  chord  C4%-blockage]  NACA 
0012  airfoil  model. 

Iteration  experiments  with  the  4%-blockage  model  at  M«,  =  0.85  and 
ot  =  1"  were  performed.  They  were  inconclusive  because  of  flow-field  unsteadi¬ 
ness.  Wall  control  was  used  to  obtain  a  first  iterative  step  toward  unconfined 
flow,  but  the  shock  wave  on  the  lower  surface  fluctuated  over  about  15%  of 
the  chord.  Subsequent  attempts  to  iterate  at  this  test  condition  did  not 
lead  to  a  steady  flow  field,  and  it  was  concluded  that  this  test  condition 
was  not  suitable  for  iteration  at  this  stage  of  the  tunnel  development.  Next, 
experiments  were  performed  with  the  4% -blockage  model  at  »  0.8  and  QU  =  4*, 
for  which  there  is  no  shock  wave  on  the  lower  airfoil  surface  and  the  upper 
shock  should  extend  almost  to  the  static  pipe.  Iterations  in  this  case  ex¬ 
hibited  steady  flow  at  each  step,  but  a  convergence  anomaly  occurred.  That 
is,  the  third  iterative  step  had  almost  converged,  but  the  fourth  appeared  to 
diverge.  In  all  of  these  experiments,  the  lack  of  a  sufficient  number  of  re¬ 
liable  normal  velocity  measurements  at  the  control  surfaces  hindered  the 
iterative  process. 

The  objectives  of  the  present  investigation,  then,  were  to  continue 
the  assessment  of  the  Calspan  two-dimensional,  perforated-wall,  segmented -plenum 
implementation  of  the  adaptive-wall  concept  at  supercritical-flow  conditions 
for  which  model- induced  shock  waves  reach  the  walls.  An  important  aspect  of 
the  investigation  was  the  development  of  a  new  instrumentation  technique  for 
determining  the  normal  velocity  component.  This  technique  involves  measure¬ 
ment  of  the  static  pressure  at  diametrically  opposed  orifices  on  the  top  and 
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bottom  of  a  static  pipe  and  relating  the  resulting  distributions  of  the  pres¬ 
sure  differences  to  the  streamwise  derivative  of  the  normal  velocity.  Details 
of  the  flow-measuring  characteristics  of  the  modified  instrumentation  and 
the  implications  on  tunnel  operation  and  data  reduction  procedures  are  described 
in  Section  2.  The  experimental  results  are  given  in  Section  3  beginning  with 
calibrations  of  the  instrumentation  without  the  model  present  in  the  tunnel. 

Next,  experiments  with  the  4%-blockage  model  at  M«,  =  O.S  andoc-  =  4*  are 
reported,  followed  by  the  results  of  experiments  with  the  same  model  at  =0.9 
andot  =  4*’.  In  the  latter  case,  the  shock  waves  extend  from  the  model  to  both 
walls.  The  investigation  is  summarized  by  concluding  remarks  in  Section  4. 

In  the  Appendix,  additional  aspects  of  static-pipe  technique  are  presented, 
namely  an  alternative  derivation  of  the  static -pipe  characteristics  to  verify 
the  completeness  of  the  analysis,  and  some  possible  alternative  implementations 
in  both  two-  and  three-dimensional  flow  fields. 
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2 . 0  INSTRUMENTATION 


2.1  BACKGROUND 

3  4  6 

In  the  experiments  prior  to  this  investigation,  '  '  local  flow  angle 
was  measured  b/  aerodynamic  probes  and  the  local  static  pressure  by  static  pipes. 
These  measurements  were  used,  respectively,  to  determine  the  normal,  V,  and 
streamwise,  u-  ,  disturbance  velocity  components  directly.  However,  the  flow- 
angle  measurements  always  limited  the  efforts.  The  inherent  limitations  in  the 
probe  technique  arise  because  of  the  small  scale,  since  blockage  considerations 
limit  the  number  of  probes  that  can  be  used.  Therefore,  it  is  difficult  to 
obtain  sufficient  measurements  to  define  adequately  the  V  distributions.  More¬ 
over,  the  small  probes  that  must  be  used  are  very  sensitive  to  contamination 
from  oil  present  in  the  air  stream  of  the  tunnel.  An  alternative  technique  is 
to  measure  the  static  pressures  at  one  control  surface  and  the  difference  between 
those  pressures  and  the  pressures  at  a  second  surface  slightly  farther  away  from 
the  model.  In  effect,  this  can  be  regarded  as  measuring  the  local  static  pressure 
and  its  gradient,  from  which  the  streamwise  derivative  of  the  normal  velocity 
can  be  inferred.  The  advantage  in  this  measurement  technique  is  that  static 

pressure  is  easy  to  measure  with  good  precision  and  one  can  easily  obtain  good 

12 

spatial  resolution.  The  TSFOIL  computer  code  was  used  to  calculate  preliminary 
estimates  of  the  pressure  differences  to  be  expected,  and  they  are  readily 
measurable  in  the  immediate  vicinity  of  the  airfoil, 


12.  ^rman,  E.M. ,  Bailey,  F.R.  and  Johnson,  M.L.  "TSFOIL-A  Con^juter  Code 

Transonic  Calculations,  Including  Wind-Tunnel 
Wall  Effects  and  Wave-Drag  Evaluation"  Paper  No.  26  in  Aerodynamic 
Analyses  Requiring  Advanced  Computers,  NASA-SP- 347-Part  2,  March  1975 
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2.2  CONFIGURATION  OF  NEW  INSTRUMENTATION 


Several  candidate  static-pipe  configurations  for  obtaining  these 
differential  measurements  were  examined.  On  the  basis  of  available  materials 
and  ease  of  construction,  it  was  decided  that  the  new  pipes  would  have  a  5/8-inch 
OD  circular  cross  section  with  diametrically  opposed  orifices  on  the  top  and 
bottom  of  the  pipes  in  the  vicinity  of  the  model.  There  are  18  pairs  of 
these  orifices  extending  9  inches  upstream  and  11  inches  downstream  of  the 
junction  between  plenum  chambers  6  and  7,  as  shown  in  Figure  1.  These  dual 
orifices  span  the  region  where  the  static  pressure  differences  are  equal  to 
or  greater  than  the  resolution  capability .  Upstream  and  downstream  of  this 
region  the  static  pipes  have  orifices  which  extend  the  full  length  of  the  test 
section  along  the  side  of  the  pipe  facing  the  model.  The  most  forward  static 
pressure  orifice  on  each  pipe  is  connected  to  a  manifold  and  the  reading  is 
taken  to  be  the  free- stream  static  pressure,  .  All  remaining  33  pressures 
on  the  model  sides  of  the  pipes  are  measured  relative  to  .  In  addition, 
the  differential  pressures  between  the  18  opposing  pairs  of  orifices  on  each 
pipe  are  measured.  The  differential  pressure  transducers  used  have  a  probable 
error  of  0.001  psi  or  less,  and  the  read-out  system  has  a  resolution  of  0.001  psi. 

After  fabrication,  the  new  static  pipes  were  mounted  with  their 
centerlines  nominally  four  inches  from  the  test  section  centerline,  as  the 
original  pipes  were.  However,  in  order  to  use  the  probe  flow-angle  data  in 
conjunction  with  the  differential  pressure  data  (as  will  be  discussed  in  the 
next  section),  the  probes  were  relocated  to  the  plane  of  the  static  pipe 
centerlines  from  their  original  locations  four  and  one-half  Inches  from  the 
test  section  centerline.  In  the  original  configuration,  four  of  the  probes 
in  the  vicinity  of  the  model  were  mounted  through  the  test  windows  and  their 
frames.  It  was  not  feasible  to  relocate  or  remove  these  probes  so  they  were 
retracted  to  lie  against  the  side  walls.  However,  two  new  flow-angle  probes 
were  fabricated  and  mounted  in  the  vicinity  of  the  model.  The  resulting 
locations  of  the  flow-angle  probes  at  the  lower  control  surface  are  shown  in 
the  side  view  of  Figure  1,  while  the  upper  probes  are  shown  in  the  top  plan- 
form  view.  The  lateral  staggering  of  the  probes  was  chosen  to  avoid  inter¬ 
ference  effects. 
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2.3  STATIC-PIPE  ANALYSIS 


The  new  static  pipes  are  used  to  determine  both  disturbance  velocity 
conqjonents.  The  static  pressures  on  the  top  and  bottom  of  each  pipe  at  a 
given  streamwise  station  are  measured  and  averaged  to  infer  the  streamwise 
disturbance  velocity  component  at  the  control  surface.  In  addition,  the  dif¬ 
ference  in  pressure  across  the  pipe  can  be  used  to  infer  the  streamwise 
gradient  of  the  normal  velocity  component  at  the  control  surface.  In  the 
following  paragraphs  the  interaction  between  the  flow  to  be  measured  and  the 
pipe  is  analyzed. 


A  schematic  of  the  model  and  the  upper  pipe  is  shown  in  Figure  2. 

The  Xrt^  coordinate  system  in  Figure  2  is  used  in  the  experiments  as  well  as 
in  the  analysis.  The  origin  is  at  the  tunnel  centerline  and  the  ^  axis  co¬ 
incides  with  the  vertical  line  passing  through  the  junctions  between  plenum 
chambers  6  and  7  and  plenum  chambers  16  and  17,  see  Figure  1.  With  the  models 
installed,  the  origin  is  0. 250'C  behind  the  leading  edge  of  the  6%-blockage 
model  and  0.175c  behind  the  leading  edge  of  the  4%-blockage  model  when  the 
models  are  at  zero  angle  of  attack. 


The  nose  of  the  pipe  is  located  well  upstream  of  the  test  section 
in  the  contraction  region.  Therefore,  the  pipe  can  be  considered  to  be  doubly 
infinite  in  length  and  of  constant  radius,  R  ,  so  far  as  the  flow  analysis  is 
concerned.  The  flow  in  the  vicinity  of  the  pipe  is  assumed  to  be  inviscid 
and  irrotational  and  to  consist  of  perturbations  about  a  free  stream  of  speed 
U„.  Accordingly,  a  perturbation  velocity  potential,  ,  exists  and  will  be 
written  as 

where  ^ the  two-dimensional  potential  of  the  disturbances  introduced 
by  the  model,  including  its  interaction  with  the  wall,  and  )  is 

the  three-dimensional  potential  of  the  disturbances  arising  from  the  inter¬ 
action  of  the  pipe  with  the  model/wall  flow  field.  The  objective  of  the 
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analysis  is  to  determine  4  and  its  derivatives  on  the  pipe  in  terms  of  the 
model /wall-induced  flow  field  that  is  to  be  determined,  namely  <p^  and  its 
derivatives.  The  disturbances  introduced  by  the  model  and  wall  are  assumed 
to  be  of  order  of  magnitude  6  ,  which  represents  the  thickness  ratio,  2*  , 
or  angle  of  attack,  o6  ,  of  the  model.  The  characteristic  length  of  the  model/ 
wall  flow  field  is  the  airfoil  chord,  c  ,  and  therefore  a  small  parameter 
representing  the  pipe  flow  field  is  .  Both  £  and  6  are  of  the  same 

order  of  magnitude  so  the  analysis  must  be  carried  out  to  include  terms  of 
order  £  ^ ,  5  *  and  £  <5“ . 

The  doubly-infinite  static  pipe  is  as  shown  in  Figure  3.  The  coordi¬ 
nates  are  the  same  as  in  Figure  2,  except  for  a  transformation  to  the  pipe 
centerline, 

f  ■  3 

The  model/wall-induced  velocity  components  along  the  pipe  centerlines,  and 
,  are  given  by 

Tr,w  * 

The  superscript  m  has  been  omitted  from  and  V,  to  simplify  the  notation, 
since  these  are  the  quantities  to  be  deduced  from  the  pipe  measurements. 
Although  the  normal  component  1/^  varies  with  over  the  pipe,  this  effect 
is  of  higher  order  and  &o  i  can  be  considered  as  the  incident  \'elocity 

normal  to  the  pipe,  as  shown  in  Figure  3.  Since  the  pipe  is  assumed  doubly 
infinite  in  length,  the  streamwise  perturbation  t^'^does  not  interact  with  the 
pipe  although  its  variation  with  c^'  is  taken  into  account  as  described  below. 
Since  typical  maximum  values  of  correspond  to  an  angle  of  only  a  few 

degrees,  it  is  satisfactory  to  use  Munk-Jones  theory  for  the  attached  flow 
over  slender  bodies  at  an  angle  of  attack, The  analysis  given  below  is  cast 
in  terms  of  conventional  s lender-body  theory.  However,  a  more  rigorous 

13.  Sears,  W.R.  "Small  Perturbation  Theory",  Section  C  of  General  Theory 
of  High  Speed  Aerodynamics,  Vol.  VI  of  High  Speed  Aerodynamics  and 
Jet  Propulsion.  Princeton  U.  Press,  1954. 


(3> 

C4) 
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analysis  using  matched  asymptotic  expansions  is  presented  in  Appendix  A.l. 

This  latter  analysis  verifies  that  the  results  obtained  are  correct  to 
second  order. 

Application  of  Munk- Jones  theory  to  this  compressible- flow  configura¬ 
tion  gives  the  perturbation  potential  for  the  pipe  as  siiaply  that  for  a 
doublet  in  incompressible  crossflow  of  strength  (X'i,  namely 

where  the  cylindrical  coordinates  of  Figure  3  are  used  for  convenience.  This 
is  the  potential  necessary  to  insure  zero  normal  flow  at  the  pipe  surface 
for  the  incident  crossflow  .  The  velocity  components  corresponding 

to  an  arbitrary  field  point,  again  in  cylindrical  coordinates,  are 


.  (6] 

Ct: ,  r,  io)  »  {!) 

vj*  ^  >  *^1  ?  =  viix)  R^  (8) 

where  (■X'i  =  cl  /cL%  .  When  these  are  resolved  into  the  X,  u  ,  ^ 

coordinates,  they  become  ® 

P 

(JX ,f',  Ui')  =  ~vj(x^  R*  Oj/r  (9) 

«  -  v^(X}  R*  2ui/r‘ 


,  r,  u; )  =  -  Cx;j 26c7/r® 


When  these  expressions  are  evaluated  on  the  pipe,  r»R  ,and  are 
added  to  the  velocity  components  of  the  mode 1/wall- induced  flow  field. 
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the  total  perturbation  velocities  u  ,  v  .cj  on  the  pipe  are  obtained.  Treat¬ 
ing  the  streaxiwise  component  first,  it  follows  that 

UL  «  tjS  (.X ,  R  t  to) 

On  the  pipe  ^ varies  with  and  the  first  two  terms  of  a  Taylor's  series  can 
be  written  as 

=.  a.<t>  a' 

where  the  terms  neglected  are  of  order  f 5*  .  If  the  equation 

*  -  /?  uj  (14) 

which  holds  on  the  pipe,  is  used  with  the  condition  for  irrotational  flow 
iSuT"/ di^*  ]  =  V/CX)  (15) 

it  follows  that 

a  ^  a  iX)  -  Vg' LX)  ^  (16) 


When  this  is  added  to  Eq.  (9),  the  result  on  the  pipe  is 
U.  «  txJX")  —  2  if’^LX')  R  CO 

O  ^ 


(17) 


where  the  second  term  has  equal  contributions  from  the  ^  -derivative  of  the 
model/wall-induced  flow  at  the  pipe  and  from  the  crossflow-induced  flow  over 
the  pipe.  Since  LX)  is  of  order  £,  and  R  is  of  order  B6.  the 

terms  in  Eq.  (17)  are  of  first  order  and  second  order,  respectively.  The 
components  V  and  ut  on  the  pipe  follow  directly  from  Eqs.  (10)  and  (11)  and 
the  fact  that  for  two-dimensional  flows  OJ^LX)^0  ,  namely 


V  =  ti;  -  tr  C5C)  2  a) 


(18) 


or  =  ~Vg(X)  2  to 


(19) 


The  exact  pressure  coefficient  for  isentropic  flow,  when  expanded 
to  second  order  in  the  perturbation  velocity  components,  is 
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t  z 

where  yS  =  f~^te  ■  It  is  convenient  at  this  point  tonondimensionalize  ac¬ 
cording  toC'^)  =  C  )/  U*  and  X  *  Z/c  .  Then,  substitution  of  Eqs.  (17} 
to  (19)  into  Eq.  (20}  and  retention  of  all  second-order  terms  in  £  and  <S  give 
for  the  pressure  coefficient  on  the  pipe, 

CpCTC^  (*i)  -  -  2  a,  <£>  +  4<S  Z-Zcz")  c<Kj  CO  -/3^ 

—  4  V  ^(X)  U> 

Therefore,  the  objective  of  the  analysis  has  been  achieved  by  relating  the 
pressure  coefficient  measured  on  the  pipe  to  the  mode 1/wall -induced  velocity 
components  a,  and  at  the  pipe  centerline. 


3  4  6 

The  static  pipes  used  in  the  earlier  experiments  '  ’  had  a  single  row 


of  orifices  at  u>=  n/z  ,  so 


Cp^  S  Cp(x,Tt/Z')  =  -  ~  -  4  vZiX) 


(22} 


In  the  earlier  experiments,  this  expression  was  approximated  by  the  leading, 
first-order  term  in  order  to  evaluate  u  from  the  Cp  measurements. 

This  linearization  was  consistent  with  the  external-flow  calculations  for 
unconfined  flow  since  they  are  based  upon  transonic  small -disturbance  theory. 


The  new  static  pipes  have  two  rows  of  orifices  located  at  a>=o  and 
TTii^- hj  and  in  Figure  3},  for  which  Eq.  (21}  gives 


C 


CpCX,o)  ~  -  3  a,  (£?  +  4(5“  'CX)  -  /3^  a/cx  J 
CpC%,7r)=  -3a^C*)-4(S 


(23} 

(24} 


If  these  are  added  and  rearranged,  it  follows  that 


a.  cx)  + 

o 


/3“a,®(;d) 


Cp,  +  Cg 


{25} 
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which  can  be  solved  to  obtain 

which  reduces  to  the  linear  relationship 


4 


(26) 


(27) 


for  small  disturbances.  If  Eqs,  (23)  and  (24)  are  subtracted  and  rearranged, 
it  follows  that 


v/c£) 


8S 


(28) 


The  result  emphasizes  the  importance  of  accounting  properly  for  the  inter¬ 
action  between  the  pipe  and  the  model/wall  flow  field.  If  that  interaction 
were  ignored  and  only  the  variation  of  with  considered,  the  factor  of 
8  in  the  denominator  would  be  reduced  to  4,  so  a  totally  incorrect  result 
would  have  been  obtained. 


The  CX)  distribution  over  the  interval  of  the  static  pipe  dif¬ 
ferential  measurements  is  found  by  integration,  namely 

£ 

f  lv/(4')loL4  t29) 

where  is  given  by  Eq.  (28)  and  V,  c£a>  is  the  value  of  at  some 

reference  point  .  The  probe  measurements  are  used  to  determine  ir^  (X^)  , 
as  well  as  to  complete  the  determination  of  the  distribution  upstream  and 
downstream  of  the  limits  of  the  Cp  -Cp^  measurements.  This  combination  of 
probe  and  pressure  gradient  measurements  provides  the  distributions  which 
are  used  as  the  boundary  conditions  for  the  external-flow  calculations. 

The  analysis  given  here  for  the  flow  about  the  static  pipe  is  no 
longer  valid  in  the  immediate  region  where  a  shock  wave  from  the  model  or  wall 
impinges  on  the  pipe.  Accordingly,  special  treatment  is  necessary  when  this 
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occurs.  The  procedure  for  handling  the  data  near  shocks  is  described  in 
Section  2.4,  where  results  for  a  typical  case  are  presented  and  discussed. 

2.4  NOilMAL-VELOCITY  DATA  PROCEDURE 

In  this  section,  the  notation  will  be  simplified  from  that  of  the 
analysis  for  use  in  the  remainder  of  the  report.  In  particular,  the  subscript 
zero  will  be  dropped  from  the  disturbance  velocity  components  u,  and  at 
the  control  surfaces  and  from  the  control  surface  location  .  Thus  u.  and 

V  will  refer  to  these  components  evaluated  at  the  control  surfaces  located 
at  -  ±\r,  . 

The  data  procedure  is  best  described  by  considering  an  example.  The 
example  chosen  is  the  second  step  in  the  experimental  iteration  for  M^,  =0.8 
and  a  =  4®  that  is  discussed  in  Section  3.2.  The  data  for  cLiV’/Ug^) 
as  calculated  by  Eq.  (28]  from  the  static-pipe  measurements,  are  presented 
in  Figures  4a  and  4b  for  the  upper  and  lower  control  surfaces,  respectively. 

The  estimated  unconfined-f low  distributions  calculated  by  TSFOIL  are  also  given 
in  Figure  4,  but  are  for  reference  only  since  there  is  no  reason  to  believe 
that  the  data  at  this  step  of  the  iteration  should  agree  with  the  estimate. 

The  data  procedure  at  the  lower  control  surface  will  be  discussed 
first  since  there  was  no  shock  wave  below  the  airfoil.  The  dL  (v/ Ug,)/U(X/c) 
data  of  Figure  4b  were  fit  by  a  cubic  spline  over  the  entire  interval  from 
7C/c=  -2.25  to  z/c*  2.75  for  which  measurements  were  made.  This  spline 
fit  was  then  integrated  according  to  Eq.  (29)  with  the  reference  velocity 
V(Z.^/C)/0^  temporarily  set  equal  to  zero  at  z/c  =  -2,25  to  give  the  dashed 
curve  shown  in  Figure  5b.  The  actual  reference  velocity  V  (X^/c)/u^v3iS 
then  established  by  a  least-squares  fit  of  the  integrated  curve  to  the  four 
probe  measurements  of  Tr/U„  at  Z/c  -  -1.48,  -1.12,  0.83  and  2. 52  in  Figure  Sb. 
The  resulting  curve  was  then  faired  upstream  through  the  first  three  probe 
measurements  and  downstream  through  the  last  two  to  arrive  at  the  final 
distribution  shown  in  Figure  Sb. 
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At  the  upper  control  surface  the  shock  wave  appears  just  to  reach 
the  upper  pipe,  as  indicated  by  Schlieren  observations.  As  shown  in  Figure  4a, 
there  are  very  large  gradients  of  ol  Cas’/O  in  the  vicinity  of  the  shock. 

Unfortunately  the  spacing  of  the  measurements  is  still  large  enough  that  the 
curve  cannot  be  clearly  resolved  in  this  region.  Therefore  separate  treat¬ 
ment  is  required  upstream  and  downstream  of  the  shock,  whose  location  in  any 
experiment  is  apparent  from  the  cL(v-/Uta)/d(x/c)^tdi  and  the  Schlieren  observa¬ 
tion.  Accordingly,  the  cL(lr/0„'i M(X/c)  data  were  fit  by  a  cubic  spline  from 
the  first  measured  point  at  X/c  =  -2.2S  to  the  last  point  before  the  shock  at 
“X/c  =  0.25.  The  spline  fit  was  then  integrated  up  to  ?^/c  ~  0.25  with  the 
reference  velocity  tr (aJ, /c)/U„  temporarily  set  equal  to  zero  at  x/c  ~  -2.25 
to  give  the  dashed  curve  in  Figure  Sa.  The  actual  reference  velocity  was 
chosen  by  adjusting  the  integrated  curve  to  pass  through  the  probe  measure¬ 
ment  at  X/c  ■  -1.95  and  the  curve  was  faired  upstream  through  the  first  two 
probe  measurements.  Downstream  of  the  shock  a  similar  procedure  was  used. 

The  c£  ('V/ ) /e^ f aJ/C  )  data  from  X/c  -  0.50  to  %/c  ~  2.75  were  spline  fit  and 
then  integrated  upstream  from  a  temporary  reference  value  of  zero  at  x/c.  ~  2,75 
to  Z/c.  »  0.50  to  give  the  other  branch  of  the  dashed  curve  in  Figure  5a,  This 
branch  was  adjusted  by  a  least- squares  fit  to  the  two  probe  measurements  at 
X/c.  »  1.25  and  2.05.  It  was  then  faired  downstream  through  the  last  three 
probe  measurements.  In  the  vicinity  of  the  shock,  the  upstream  and  downstream 
branches  were  simply  extrapolated  until  they  intersect  to  give  the  final 
distribution  in  Figure  5a. 

It  is  clear  that  additional  measurements  in  the  vicinity  of  the 
shock  wave  would  be  highly  desirable,  but  this  was  not. possible  with  the 
present  configuration.  In  fact,  it  would  be  difficult  to  construct  a  pipe 
with  sufficient  measurement  orifices  to  cover  all  shock  locations  for  a  wide 
range  of  models  and  operating  conditions.  Accordingly,  an  alternative  such 
as  pipes  which  translate  streamwise  offers  a  better  possibility  for  more  com¬ 
plete  definition  of  the  distributions  in  future  applications.  Despite  this 
limitation  of  the  existing  configuration,  the  new  static  pipes,  in  conjunction 
with  the  probes,  provide  an  improved  technique  for  determining  the  normal- 
velocity  distributions  compared  to  that  of  the  probes  alone. 
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2.5  ADAPTIVE-WALL  OPERATING  TECHNIQUE 

The  adapt ive-ua 11  operating  technique  is  basically  the  same  as  that 
used  in  the  past,  with  data  procedures  modified  to  suit  the  new  static -pipe 
configuration.  An  initial  approximation  to  the  static  pressure  distributions 
on  the  pipes,  say,  is  found  by  using  TSFOIL  to  estimate  the  unconfined- 

flow  distributions  of  0-/(4  and  Cc-'cJ  for  use  in  evaluating  Eq.  (23). 

The  resulting  distributions  are  set  in  the  tunnel  and  the 

distributions  are  measured  along  with  the  probe  readings.  From  these  measure¬ 
ments,  the  actual  experimental  distributions  are  evaluated  as  described 

above  and  used  as  boundary  conditions  for  the  external-flow  calculation  of 
the  corresponding  unconfined- flow  distributions  64  .  These  calculations 

are  made  using  a  computer  code  based  on  finite-difference  solutions  to  the 

4  6 

transonic  small-disturbance  equations.  This  code  had  been  developed  and  used 
in  the  earlier  investigations.  The  new  calculated  distributions 
are  used,  in  turn,  along  with  aLCv'^i^ Ufg)/d(X/c'i  and  Eq.  (23)  to  determine 
calculated  distributions,  Cp^^  say,  which  are  relaxed  to  get  the  distribu¬ 
tions  to  be  set  in  the  tunnel  in  the  second  iterative  step,  namely 

-  Cl  -  i)  Cp\  -Ik  (30) 

where  is  the  relaxation  factor.  After  setting  ,  the  appropriate 

measurements  are  repeated  and  the  iteration  continued  until  convergence  is 
achieved. 
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3.0  EXPERIMENTS 

3.1  INSTRUMLNIATION  CALIBRATIONS 

The  first  experiments  were  a  series  of  tests  without  the  airfoil 
model  installed  in  the  tunnel.  For  convenience,  these  will  be  referred  to  as 
empty- tunnel  tests.  These  tests  were  necessary  in  order  to  calibrate  the  two 
new  flow-angle  probes  over  the  geometric  angle-of-attack  range  from  -4"  to 
>4'*  at  free-stream  Mach  numbers  from  0.5S  to  0.95.  It  was  found  that  the 
response  of  the  new  probes  to  changes  in  angle  of  attack  is  linear,  which  is 
in  agreement  with  the  response  of  the  other  flow- angle  probes.  These  experi¬ 
ments  also  were  required  to  provide  new  probe  readings  at  zero  flow  angle  over 
the  same  Mach  number  range,  since  all  the  probes  had  been  relocated  in  the 
plane  of  the  static-pipe  centerlines.  Later,  during  the  course  of  the  itera¬ 
tion  experiments,  additional  empty-tunnel  runs  were  made  to  acquire  updated 
flow-angle  probe  calibrations  at  zero  angle  of  attack. 

As  a  by-product  of  these  empty-tunnel  calibrations  runs,  the  quality 
of  the  uniform  flow  in  the  test  section  was  measured.  Typical  results  at 
M^  =  0.85  are  shown  in  Figure  6  and  compared  with  earlier  data  obtained  with 
the  original  1/2-inch  OD  static  pipes.  The  results  are  presented  as  the  stream- 
wise  disturbance  velocity  component  plotted  as  a  fimction  of  the  distance 
from  the  beginning  of  the  controlled  part  of  the  test  section.  The  quality  of 
the  uniform  flow  obtained  in  the  recent  tests  matches  that  of  the  earlier 
results.  Any  variations  in  velocity  exceeding  0.5%  occur  either  near  the  begin¬ 
ning  of  the  test  section,  or  well  downstream  near  the  end  of  the  test  section. 
For  the  data  shown  in  Figure  6,  the  mean  values  of  are  -0.0008  and  -0.0009 

for  the  present  and  previous  tests,  respectively.  The  corresponding  standard 
deviations  for  both  sets  of  data  are f 0.0024. 

The  new  static  pipe  data  have  been  analyzed  for  each  of  the  test 
Mach  numbers.  The  mean  values  of  and  the  standard  deviations  are  pre¬ 
sented  in  Table  I  for  representative  cases.  The  deviation  of  the  mean  values 


23 


AEDC-TR-BO-34 


from  zero  is  significantly  less  than  the  standard  deviations  (with  one  ex¬ 
ception  -  the  upper  pipe  at  M„  =  0.55),  The  flow  quality  at  Mach  numbers  of 
O.SSj  0.75  and  0.85  is  comparable  and  the  standard  deviations  are  all  about 
-0.0020.  At  MfQ  -  0.95,  the  standard  deviation  of  all  the  data  is  approximately 
±0.0040;  however,-  this  is  due  to  deviations  at  the  far  upstream  and  downstream 
ends  of  the  test  section.  Over  the  central  20- inch  distance  near  the  model 
location,  the  standard  deviation  is  ±0.0022,  which  is  in  good  agreement  with 
the  lower  Mach  number  results. 

The  pressure  differentials  measured  by  the  opposing  pairs  of  orifices 
have  also  been  analyzed.  The  mean  values,  expressed  as  Cp~C^^  ,  and  the 
standard  deviations  are  shown  in  Table  II  for  the  same  representative  cases. 
These  orifices,  which  are  located  as  shown  in  Figure  1,  exhibit  mean  values  of 
about  ±0.0010  with  standard  deviations  of  ±0.0020.  These  mean  values  are 
small  and  comparable  to  those  for  the  overall  flow,  as  in  Table  I  and  Figure  6. 
However,  the  fact  that  the  means  are  not  zero  has  a  significant  implication 
on  the  determination  of  Tr/Ug^  with  the  pipes.  The  empty-tunnel  differential 
pressure  data  were  reduced,  spline  fit  and  integrated  according  to  Eq.  (29) 
over  the  total  length  of  the  measurements  from  2£/c  =  -2. 25  to  2.75  using  the 
standard  data  reduction  procedure.  The  resulting  cumulative  effect  of  the 
nonzero  pressure  differentials  of  one  predominant  sign  was  very  large.  This 
can  be  seen  in  Figure  7,  where  the  results  at  the  upper  pipe  for  two  different 
empty- tunnel  runs  at  M*,  =0.8  are  presented.  Since  this  plot  is  to  the  same 
scale  as  Figure  5,  a  comparison  shows  just  how  large  the  effect  is.  All 
empty-tunnel  data  were  examined  in  this  way.  It  was  found  that  the  upper 
pipe  data  integrated  in  all  cases  to  give  positive  v/U^,  distributions  of 
roughly  the  magnitude  shown  in  Figure  7.  At  the  lower  pipe,  however,  the 
results  were  different  and  contradictory.  That  is,  some  cases  integrated  to 
give  positive  "v/Ugg  (of  smaller  magnitude  than  in  Figure  7)  and  some  to  give 
negative  T//i^  .  Indeed,  at  the  lower  pipe  w'as  positive  in  one  run  and 
negative  in  the  other  for  the  experiments  shown  in  Figure  7. 
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Several  possible  causes  of  this  phenomenon  were  examined.  Calcula¬ 
tions  indicated  that  the  pipe  does,  not  deflect  statically  an  amount  sufficient 
to  cause  distributions  of  /d(z/c)Qf  the  magnitude  that  was  measured. 

This  was  verified  by  careful  measurements  along  the  length  of  the  pipes  of 
the  distance  from  the  pipe  to  the  wall.  However,  the  deflection  under  tunnel 
operating  conditions,  including  dynamic  loads  and  thermal  effects,  could  not 
be  assessed  readily.  The  normal  velocities  implied  at  both  pipes  by  the 
measurements  for  each  run  can  be  interpreted  in  terms  of  an  apparent  effective 
area  change  in  the  flow.  Such  an  area  change  implies  a  variation  in  u./1/a, 
along  the  streamwise  direction  by  consideration  of  the  continuity  equation. 

For  two  representative  cases,  the  implied  a/U,,  distribution  due  to  the  v/U„ 
measured  at  the  two  pipes  was  calculated  by  means  of  a  one-dimensional  flow 
model.  The  calculated  u/U^  values  were  very  much  larger  in  magnitude  and 
varied  with  X/c  in  an  entirely  different  way  from  the  measured  u/Uf^  values , 
e.g.,  those  shown  in  Figure  6.  Therefore  it  was  concluded  that  the  flow  is  in 
fact  uniform,  as  indicated  by  the  small  measured  u/U^a  variations,  and  the 
anomalies  in  the  pipe  differential  measurements  are  an  artifact  of  the  detailed 
geometry  of  the  individual  pipes,  since  the  characteristics  of  each  differ. 
Consequently,  this  effect  at  the  upper  pipe  was  treated  as  a  tare  reading 
to  be  Subtracted  at  each  orifice  from  the  measured  d(-1r/U„>/diX/C). 

The  integrated  v/Uf^  from  the  assumed  tare  is  shown  in  Figure  7  and  appears 
to  provide  a  reasonable  approximation  to  the  experimental  curves  for  both 
empty-tunnel  runs  at  this  .  No  tare  was  used  at  the  lower  pipe  because 
it  averaged  nearly  to  zero  for  all  runs  considered. 

3,2  EXPERIMENTS  AT  -  0.8  AND  a  =  4® 

These  experiments  with  the  4%-b lockage  model  were  performed  in 
order  to  investigate  further  the  convergence  anomaly  that  had  been  observed 
in  the  earlier  experiments,  as  mentioned  in  the  Introduction.  Initially,  this 
case  was  chosen  because  it  provides  a  flow  field  in  which  only  the  upper  air¬ 
foil  surface  is  supercritical  with  the  shock  extending  to  the  vicinity  of 
the  static  pipe,  but  not  to  the  wall.  Rather  than  beginning  where  the  earlier 
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experiments  ended,  it  was  more  convenient  with  the  new  instrumentation  to 
begin  the  iteration  over.  Accordingly  the  procedure  outlined  in  Section  2.5 
was  used. 


Some  difficulty  was  encountered  in  eliminating  an  apparent  cross- 
flow  condition  upstream  of  the  model  at  the  first  three  plenum  chambers  on 
the  upper  and  lower  walls.  This  same  difficulty  had  occurred  in  the  divergent 
fourth  iterative  step  of  the  earlier  investigation.  It  was  found  that  this 
tendency  could  be  avoided,  or  minimized,  if  the  valve  settings  at  the  up¬ 
stream  plenum  chambers  were  set  initially  at  the  empty-tunnel  settings  for 
the  same  free- stream  Mach  number.  Accordingly  a  reasonable  first  step  in  the 
iteration  was  set.  The  resulting  pressure  distribution  on  the  airfoil  was 

C 

in  only  fair  agreement  with  the  Eight-Foot  Tunnel  data,  see  Figure  8,  with 
the  shock  too  far  back.  It  should  be  recalled,  however,  that  the  Eight-Foot 
Tunnel  data  were  obtained  at  a  Reynolds  number  of  1.0  x  10^  based  on  model 
chord  while  the  present  data  with  the  4%-blockage  model  correspond  to  a 
Reynolds  number  of  0.67  x  10^.  Earlier  experiments  had  shown  in  a  short  dura¬ 
tion  run  that  at  =0.8  and  oc,-  4“,  the  shock  wave  moved  forward  as  the 
Reynolds  number  was  increased  from  0.67  x  10^  to  1.0  x  10^. 

Although  there  seemed  to  be  a  slight  residual  negative  crossflow 
upstream  from  the  probe  measurements,  nevertheless  the  pipe  differential -pres¬ 
sure  and  probe  data  were  reduced  as  described  in  Section  2.4  to  obtain  the 

distributions  to  be  used  as  boundary  conditions  for  the  external  flow 
calculations.  The  resulting  unconfined-flow  distributions  corresponding  to 

are  compared  with  the  measured  distributions  in 
Figure  9.  There  are  significant  differences  between  the  distributions  and 
they  arise  in  part  from  the  crossflow.  Following  the  adaptive-wall  iterative 
technique  of  Section  2.5,  the  results  corresponding  to  the  computed 

values  of  Figure  9  were  relaxed  according  to  Eq.  [30)  with  a  relaxation 
factor  of  0.1.  That  is,  only  very  small  changes  from  the  initial  settings 
were  called  for. 
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The  crossflow  disappeared  while  setting  the  second  step  in  the 
iterative  procedure.  This  can  be  seen  from  the  distributions  that  were 

described  in  Section  2.4  and  are  shown  in  Figure  5.  The  pressure  distribution 
on  the  airfoil,  Figure  10,  is  greatly  improved  over  that  of  the  first  itera¬ 
tive  step  in  Figure  8,  although  the  difference  in  Reynolds  number  still 
exists.  Evaluation  of  the  external-flow  calculations  using  the  ii;;„/CjBdis- 
tributions  of  Figure  5  yields  the  a  distributions  given  in  Figure  11. 

The  agreement  is  much  better  than  in  the  first  iterative  step,  especially 
near  the  model  but  there  are  still  appreciable  differences  upstream  and  down¬ 
stream. 


It  was  intended  to  iterate  further  by  taking  a  third  step,  also  with 
■^  =  0.1.  At  this  time,  however,  an  accident  occurred  while  running  the  One- 
Foot  Ttmnel  and  severe  overheating  resulted.  An  upstream  segment  of  the 
adaptive-wall  test  section  was  damaged  and  required  rebuilding.  Thus,  the 
model  and  instrumentation  had  to  be  removed.  After  the  repairs  were  effected 
and  additional  tunnel  calibrations  were  carried  out,  the  iterations  could  not 
be  resumed  easily  from  the  point  where  they  had  ended.  Nevertheless,  the 
results  of  these  two  iterations  established  the  new  measurement  technique  and 
also  provided  a  basis  of  valve  settings  for  experiments  at  higher  Mach  numbers 
where  both  control  surfaces  are  supercritical. 

3.3  EXPERIMENTS  AT  =  0.9  AND  0£  =4" 

These  experiments  with  the  4%-b lockage  model  were  performed  to 
examine  the  control  available  when  both  walls  are  supercritical.  Existing 
numerical  results  for  the  flow  fields  of  the  NACA  0012  section,  as  calculated 
by  TSFOIL,  were  analyzed  to  select  suitable  cases.  An  additional  constraint 
on  the  selection  was  the  desire  to  choose  an  angle  of  attack  sufficiently 
large  to  avoid  the  very  small  lift-curve  slope  that  is  characteristic  of 
this  airfoil  at  =  0.85  and  0.9  for  angles  of  attack  of  three  degrees  or 
less,  see  Figure  5  of  Reference  5  or  Figure  3  of  Reference  6.  This  apparently 
results  from  the  shock-wave  locations  on  the  upper  and  lower  control  surfaces. 
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As  had  been  observed  in  Reference  6,  unsteady  shock  locations  existed  when  the 
4%-blackage  model  was  tested  in  this  range,  namely  at  =  0.85  and  cc=  l”. 
Although  it  is  not  known  if  th'e  small  lift-curve  slope  and  the  unsteadiness 
are  related,  it  was  desired  to  avoid  this  angle-of -attack  range  for  the  next 
set  of  experiments.  Selection  of  M^o  =  0.9  and  a  =  4®  satisfied  the  criteria 
for  supercritical  walls  and  a  reasonable  lift-curve  slope. 

In  the  TSFOIL  calculations,  the  lift  coefficient  measured  in  the 
Eight-Foot  Tunnel  tests  was  specified  instead  of  the  Kutta  condition  at  the 
trailing  edge.  This  was  done  so  that  the  peaks  of  the  computed  distribu¬ 
tions  and  the  computed  distributions  far  upstream  and  downstream  of  the  model 
would  be  representative  of  the  actual  experimental  conditions  for  unconfined 
flow.  The  drawback  of  specifying  the  lift  in  this  Mach  number  range  is  that 
there  can  be  an  appreciable  flow  around  the  trailing  edge,  which  gives  rise 
to  obvious  errors  in  the  pressure  distribution  on  the  model-  These  errors 
are  not  of  major  concern,  but  in  some  examples,  their  effect  does  appear 
in  the  induced  velocities  at  the  control  surfaces  above  and  below  the  trailing 
edge.  Therefore,  the  calculated  TSFOIL  distributions  for  were  used 

only  as  a  guide  for  setting  the  flow  in  the  first  iterative  step  and  it  was 
not  required  that  they  be  set  exactly.  In  subsequent  iterative  steps,  of 
course,  the  measured  distributions  and  the  distributions  from  the  external- 
flow  calculations  determine  the  next  settings,  according  to  Eq.  (30^. 

In  the  initial  experiments  to  set  this  flow  field,  most  of  the 
plenum-pressure  control  valves  were  opened  to  settings  representative  of 
the  experiments  at  =0.8  and  a  =  4®.  This  was  in  accordance  with  the 
procedure  developed  in  Reference  6  for  setting  supercritical  flows  by  proceeding 
upward  in  from  a  related  subcritical  case.  However,  the  valves  at  the 
first  three  plenum  chambers  at  the  upper  wall  and  the  first  two  chambers  at 
the  lower  wall  were  set  at  the  empty- tunnel,  uniform-flow  settings  for  =  0.9. 
As  mentioned  in  Section  3.2  this  procedure  had  been  reasonably  effective  in 
eliminating  crossflow  upstream. 
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When  the  tunnel  was  run  at  these  valve  settings,  considerable  dif¬ 
ficulty  was  encountered  in  reaching  the  desired  free-stream  Mach  number.  The 
required  tunnel-drive  power  increased  very  rapidly  and  limited  the  achievable 
to  about  0.89.  Schlieren  observations  indicated  that  the  shock  waves 
on  both  airfoil  surfaces  were  at  the  trailing  edge  and  were  steady.  It  was  sus¬ 
pected  that  the  tunnel  was  choked  downstream  of  the  model  and  this  was  verified 
by  the  measured  at  the  static  pipes  and  shown  in  Figure  12.  These  re¬ 

sults  are  reminiscent  of  the  first  attempt  to  set  supercritical  flows  in 
earlier  experiments,  see  Figures  S  and  6  of  Reference  6.  That  is,  the 
distributions  were  reasonable  upstream  of  the  model,  but  the  flow  was  super¬ 
sonic  downstream.  It  should  be  emphasized  that  no  valve  adjustments  had  been 
made  at  this  point. 

The  second  run  began  with  the  same  valve  settings  and  an  of  0.89. 
Valve  adjustments  were  then  carried  out.  Increasing  the  suction  at  the  upper 
plenim  chamber  just  ahead  of  the  model  Cplenum  No.  6  in  Figure  1)  and  the 
corresponding  lower  plenum  chamber  (plenum  No.  16  in  Figure  1)  served  to  un- 
choke  the  flow.  The  power  requirements  immediately  dropped  significantly 
along  with  key  tunnel -drive- system  temperatures  and  the  flow  became  noticeably 
quieter.  The  free-stream  Mach  number  was  increased  to  0.9  without  difficulty. 

The  distributions  were  measured  without  further  adjustments  and  are 

presented  in  Figure  13.  The  flow  downstream  of  the  model  was  mostly  subcritical, 
but  still  far  from  that  desired. 

A  third  run  was  made  in  which  additional  valve  adjustments  were 
made  in  the  usual  fashion.  Upstream  of  the  model,  the  distributions 

were  set  closely  to  the  TSFOIL  values,  but  just  downstream  of  the  model  they 
were  not.  It  is  in  this  neighborhood,  however,  where  the  TSFOIL  results  are 
questionable.  The  measured  distributions  are  given  in  Figure  14  together 

with  the  u  /14b  J /li,  distributions  calculated  from  the  external  flow.  The 
large  differences  between  u.„/U„ajid  a£W;^/i4]/i4,  at  the  upper  control  surface 
upstream  and  downstream  of  the  model  arise,  in  part  at  least,  from  appreciable 
scatter  in  the  probe  measurements  beyond  the  interval  where  the  static-pipe 
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differential  pressure  measurements  are  made.  The  scatter  made  interpolation 
difficult.  The  measured  pressure  distribution  on  the  airfoil  is  compared  with 
the  Eight-Foot  Tunnel  data  in  Figure  15.  The  agreement  is  generally  satis¬ 
factory  although  the  upper-surface  shock  lies  ahead  of  its  unconfined- flow 
location  and  the  pressure  rise  associated  with  the  shock  extends  over  a  greater 
portion  of  the  chord.  The  Reynolds  number  is  still  two-thirds  of  that  for 
the  Eight-Foot  Tunnel  experiments  but  its  effect  on  the  flow  has  not  been 
established  for  these  flow  conditions.  The  results  shown  in  Figure  14  provide 
the  basis  for  a  second  iterative  step  according  to  Eq.  C30),  again  with  =  0.1. 
Although  subsequent  iterative  steps  were  not  made  for  this  case,  sufficient 
wall  control  appears  to  be  available  to  continue  the  iteration  procedure. 
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4.0  CONCLUDING  REMARKS 


The  purpose  of  this  research  was  to  continue  the  assessment  of  the 
Calspan  adaptive-wall  wind  tunnel  for  two-dimensional  flows  with  supercritical 
walls.  An  important  aspect  of  the  research  was  the  investigation  of  a  static- 
pipe  measuring  technique  for  determining  the  normal  velocity  distributions. 

This  technique  was  developed  to  overcome  the  limitations  in  the  number  of 
flow-angle  probe  measurements  that  could  be  made  in  previous  experiments. 

The  flow  about  static  pipes  in  the  presence  of  model/wall-induced 
disturbances  was  analyzed  within  the  framework  of  conventional  s lender-body 
theory  and  the  theory  of  matched  asymptotic  expansions.  This  analysis  relates 
differential  pressure  measurements  across  the  pipe  to  the  streamwise  deriva¬ 
tive  of  the  normal  velocity  at  the  pipe  centerline.  The  derivative  is  integrated 
to  obtain  the  distribution  of  the  normal  velocity  if  at  least  one  independent 
normal  velocity  measurement  is  made.  In  addition,  the  average  of  the  pressure 
measurements  across  the  pipe  determines  the  streamwise  velocity  component  at 
the  pipe  centerline.  Two  suitable  static  pipes  were  constructed.  For  uniform, 
parallel  flow  in  the  empty  test  section,  the  differential  pressure  measurements 
indicated  an  apparent  gradient  in  normal  velocity  at  the  upper  pipe  for  all 
Mach  numbers  tested,  it  has  been  concluded  that  this  gradient  is  an  artifact  of 
the  pipe  construction  and  installation.  Accordingly,  it  has  been  eliminated 
by  treating  the  empty-tunnel  gradient  as  a  tare  reading.  The  new  static-pipe 
technique,  in  conjunction  with  flow-angle  probes,  was  applied  to  cases  with 
a  model  present  at  =0.8  and  «  =  4".  Results  indicated  that  the  procedure 
is  an  improvement  over  probes  alone  and  it  was  used  in  the  subsequent  experi¬ 
ments.  Also,  static  pipes  can  be  used  to  determine  the  velocity  normal  to  the 
pipe  in  three-dimensional  model/wall-induced  disturbance  fields.  However, 
further  development  to  provide  additional  measurements  in  regions  of  large 
normal- velocity  gradients  should  be  carried  out  concurrently  with  analysis  of 
the  interaction  between  a  shock  wave  and  the  static  pipe. 
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A  series  of  tunnel-empty  runs  in  uniform,  parallel  flow  served  to 
calibrate  the  instrumentation  over  a  range  of  Mach  numbers.  In  addition, 
analysis  of  the  results  showed  that  the  flow  quality  matched  that  of  earlier 
experiments. 

The  4%-blockage  NACA  0012  airfoil  model  was  run  at  Mo,  »  0.8  and  02.  =  4*, 
and  two  converging  steps  in  an  experimental  iteration  were  accomplished.  From 
runs  made  prior  to  beginning  the  actual  iteration,  it  was  concluded  that  cross - 
flow  conditions  upstream  of  the  model  could  be  minimized  in  the  first  itera¬ 
tive  step  by  initiating  the  wall  control  adjustment  with  the  upstream  valves 
set  at  their  tunnel-empty  positions  at  that  Mcb  .  From  these  successful 
experiments,  sufficient  information  was  obtained  to  proceed  to  a  flow  with 
both  walls  supercritical.  Accordingly,  the  same  model  was  tested  at  the  same 
oe  ,  but  with  Mo  increased  to  0.9. 

Initial  experiments  at  =  0.9  and  a  =  4®  were  performed  with  the 
upstream  control  valves  at  their  tunnel-empty  settings  and  the  remaining 
valves  at  their  M*  =  0.8  and  a  »  4*  settings.  The  measured  streamwise  ve¬ 
locity  component  at  the  control  surfaces  exhibited  a  reasonable  distribution 
from  the  beginning  of  the  test  section  to  the  vicinity  of  the  airfoil  trailing 
edge  with  these  settings.  Downstream  of  this  point,  however,  the  flow  was 
choked  and  supersonic  flow  persisted  to  the  end  of  the  test  section.  The 
choking  was  relieved  by  Increasing  the  suction  at  the  upper  and  lower  plenum 
chambers  just  upstream  of  the  model.  Further  control  adjustments  led  to 
reasonable  streamwise  velocity  distributions  along  the  entire  length  of  the 
control  surfaces.  Therefore,  a  first  iterative  step  had  been  achieved,  and 
it  has  been  concluded  that  sufficient  control  is  available  to  iterate  toward  un¬ 
confined  flow  at  this  test  condition. 
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SIDE  VIEW 


Figure  1  SCHEMATIC  OF  STATIC  PIPES  AND  FLOW  ANGLE  PROBES 
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Figure  2  SCHEMATIC  OF  MODEL  AND  STATIC  PIPE 
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t  i/^Cx) 


Figures  MODEL/WALL  INDUCED  FLOW  INCIDENT  TO  THE  STATIC  PIPE 
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pO.2 


Figure  4  STATIC  PIPE  MEASUREMENTS  OF  STREAMWISE  DERIVATIVE  OF  NORMAL 
VELOCITY  COMPONENTS.  =  0.8.  a  =  4°,  4%  BLOCKAGE,  SECOND 
ITERATIVE  STEP 
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(a)  UPPER  CONTROL  SURFACE,  h/c  >  1.0 

Figure  5  NORMAL  VELOCITY  COMPONENTS  DETERMINED  FROM  STATIC  PIPES 
AND  FLOW-ANGLE  PROBES,  M„  =  0.8,  a  =  4°.  4%  BLOCKAGE,  SECOND 
ITERATIVE  STEP 
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(b)  LOWER  CONTROL  SURFACE,  h/c  =  -1.0 


Figure  5  <Cont.)  NORMAL  VELOCITY  COMPONENTS  DETERMINED  FROM  STATIC  PIPES 
AND  FLOW-ANGLE  PROBES.  M„  =  0.8,  a  =  4°  4%  BLOCKAGE,  SECOND 
ITERATIVE  STEP 
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Figure  7  APPARENT  NORMAL  VELOCITY  COMPONENT  MEASURED  BY  UPPER  STATIC 
PIPE  IN  UNIFORM,  PARALLEL  FLOW  WITHOUT  MODEL,  =  0.8.  h/c  =  1.0 
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u[v„/U^]/U„  FROM 
EXTERNAL-FLOW  CALCULATION 


«0»05 


(bl  LOWER  CONTROL  SURFACE,  h/c- -1.0 
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Figure  11  STREAMWISE  DISTURBANCE  VELOCITY  COMPONENTS,  M ^  =  0.8,  o  =  4° 
4%  BLOCKAGE,  SECOND  ITERATIVE  STEP 
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x/c 

la)  UPPER  CONTROL  SURFACE,  h/c  =  1.0 

Figure  12  MEASURED  STREAMWISE  DISTURBANCE  VELOCITY  COMPONENTS,  M  =  0.9 
a  »  40,  4%  BLOCKAGE,  INITIAL  VALVE  SETTINGS  ” 
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Figure  13  MEASU  R  ED  STREAMWISE  DISTURBANCE  VELOCITY  COMPONENTS,  M 
a  =  4®,  4%  BLOCKAGE,  PRELIMINARY  VALVE  ADJUSTMENTS 
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-  uIv„/U„]/U„  FROM 

EXTERNAL-FLOW  CALCULATION 
•  CONTROL  POINT 
O  MEASURED  DATA 


Figure  14  STREAMWISE  DISTURBANCE  VELOCITY  COMPONENTS,  M„  “0.9, 
a  =  4®,  4%  BLOCKAGE.  FIRST  ITERATIVE  STEP 
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-  u  [v„/U^I/U„  FROM 

EXTERNAL-FLOW  CALCULATION 
•  CONTROL  POINT 

O  MEASURED  DATA 


x/e 

(b)  LOWER  CONTROL  SURFACE,  h/c  -  -1.0 


Figure  14  (Cont.)  STREAMWISE  DISTURBANCE  VELOCITY  COMPONENTS,  -  0.9, 
a  =•  40  4%  BLOCKAGE,  FIRST  ITERATIVE  STEP 
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Figure  15  COMPARISON  OF  MEASURED  AIRFOIL  PRESSURE  DISTRIBUTIONS, 
M«  =  0.9,  a  =  4®,  FIRST  ITERATIVE  STEP 
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TABLE  I 

EMPTY  TUNNEL  FLOW  MEASUREMENTS  OF  M,/  U„ 


MACH 

UPPER  PIPE 

LOWER  PIPE 

BOTH 

PIPES 

NUMBER 

_  MEAN 

STD.  DEV. 

MEAN 

STD.  DEV. 

MEAN 

STD.  DEV 

.55 

-0.0017 

t 0.0016 

-0.0003 

±  0.0021 

-0.0010 

±0.0020 

.75 

-0.0008 

10.0016 

-0.0006 

10.0017 

-0.0007 

±0.0016 

.85 

-0.0008 

±0,0025 

-0.0008 

±  0.0024 

-0.0008 

±0.0024 

.95 

-0.0000 

±0.0046 

-0.0014 

±0.0036 

-0.0007 

±0.0042 

TABLE  II 

DIFFERENTIAL  MEASUREMENTS  ACROSS  PIPE  OF  -  Cp^ 
IN  EMPTY  TUNNEL 


MACH 

UPPER  PIPE 

LOWER  PIPE 

NUMBER 

MEAN 

STD.  DEV. 

MEAN 

STD.  DEV. 

.55 

+0,0004 

±0.0024 

+0.0008 

±0.0022 

.75 

+0,0014 

±0.0024 

-0.0010 

±0.0016 

.85 

+0.0014 

±0.0026 

-0.0010 

±0.0018 

.95 

*0.0016 

±0.0030 

-0.0012 

+0.0030 
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APPENDIX 


Several  features  of  the  static-pipe  analysis  and  applications  are 
presented  here.  The  first  is  a  more  rigorous  derivation  of  the  static  pipe 
analysis  using  matched  asymptotic  expansions.  This  derivation  served  to 
verify  that  the  slender-body  analysis  in  Section  2.3  is  correct  and  consistent 
to  second  order.  The  second  is  an  alternative  application  of  a  static  pipe 
in  tvro  dimensions  and  the  third  is  a  static-pipe  application  in  a  three-dimen¬ 
sional  disturbance  field. 

A.l  STATIC-PIPE  ANALYSIS  BY  MATCHED  ASYMPTOTIC  EXPANSIONS 

The  geometry  of  the  flow  is  given  in  Figures  2  and  3,  as  before. 
Again,  the  flow  is  ass\imed  to  be  subsonic,  but  compressible,  so  that  the 
Prandtl-Glauert  approximation  applies.  The  total  velocity  potential 
is  of  the  form 

^  (Al) 

where  i  =  X/c  ,  etc.  as  in  the  main  text,  and 

-  4>"'/  0^  e 
^2  ^  4>'‘/U„€ 

in  terras  of  the  main  text  notation  of  Eq.  Cl]*  where  ^"^is  the  two  dimen¬ 
sional  perturbation  potential  of  the  disturbances  introduced  by  the  model 
and  wall  and  (^^is  the  three-dimensional  perturbation  potential  arising  from 
the  interaction  of  the  pipe  with  the  model/wa 11- induced  flow  field.  The 
surface  of  the  pipe  is  described  by 

where  ,  the  nondimens ional  pipe  radius.  In  the  Prandtl-Glauert  ap¬ 

proximation,  the  potential  satisfies 
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f?  =  0  (A2) 

subject  to  the  boundary  condition 

V<p^'  ^ si  •  O  on  ^  =  C  (A3) 


As  in  the  main  text,  it  is  assumed  that  6  is  a  small  parameter,  that 
is,  6  <<  }  ,  and  the  technique  of  matched  asymptotic  expansions^  is  used  to  find 
the  appropriate  expansion  for  <P^  .  As  <5 — ►O  ,  the  pipe  shrinks  to  a  line. 
Tentatively,  then,  the  outer  solution  is  represented  by  a  line  doublet  distribu¬ 
tion  plus  a  line  source  distribution,  viz. 


r  P(XJ(X-X,)  d 

^  A 


<r,iS}  r* 

zft  J  [(x-zy*/3’rn''^ 

«  60 


CA4) 


where  ^  ^  '*'  }■  •  Th®  doublet  strength  P  will  be  determined 

by  matching  with  an  inner  solution  to  first  order  in  <5  ■  The  source  term 

(second  integral  expression)  is  not  needed  to  first  order  in  ,  so  both  its 

14 

strength  ^  Ot*)  and  gauge  function  (T^  CS)  will  be  determined  by  matching 
with  higher-order  inner  solutions. 


A. 1.1  Inner  Solution 


The  inner  problem  is  obtained  by  stretching  the  coordinate  normal 
to  the  pipe  by 

Y  -  ^'/S  ,  e  .|/f 


14.  Van  Dyke,  M.  Perturbation  Methods  in  Fluid  Mechanics,  Academic  Press, 
New  York,  ISfiT! 
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Then  must  be  expanded  accordingly  as 

(p,(^,6y)  -  +  S  <^„(x)y  +  6^<l>,g(x)y^  +  0(6")  [A53 

It  should  be  noted  that  ti  4>io^^^/dX  is  proportional  to  the  desired  streamwise 
disturbance  velocity  component  along  the  pipe  centerline  and  <S  t^^^C^lis 
proportional  to  the  normal  con^onent  there. 

Next,  a  total  inner  solution  is  assumed  of  the  form 

+  a  (A6) 

It  can  be  shovm  that  both  ^ ^  and  satisfy  the  two-dimensional  Laplace's 

equation  in  the  crossflow  plane,  while  the  boundary  condition,  Eq.  (A3), 
reduces  to 


and 


both  on 


2  +  2  2 

-  -2V'0„ 

CA7) 

+  22^3^ 

(A8) 

Y^+  2®.  7  -  O 

Physically,  it  may  be  seen  that  the  problem  for  corresponds  to 
that  of  a  cylinder  with  unit  radius  immersed  in  a  flow  of  magnitude  (p„  . 

The  corresponding  solution  is 

^2  “  <P„Y/(y%z‘)  -f-  (A9) 

or,  in  the  polar  coordinates  of  Figure  3  {  Z  “  -  (R  co  ) , 

”  *^it  m/  IR.  +  G-jd) 
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This  equation  is  of  the  same  form  as  Eq.  (55  in  the  derivation  in  the  main 
text.  Matching  with  the  outer  solution  will  be  carried  out  in  Section  A. 1.2 
to  obtain  <3-,  C^^and  Ptx.)  . 

The  solution  for  is  most  easily  found  using  polar  coordinates, 
in  terras  of  which  Eq.  (A8)  becomes 

^3^  -  -  2  4*, 2  (All) 

on  =  1,  where  the  subscript  denotes  differentiation  with  respect  to  rfZ . 
After  suitable  manipulation,  Poisson's  integral  formula  can  be  used  to  solve 
for  $3.  First,  let 

^3  ”  *  K  GgCx)  (AI25 

where  K  is  a  constant  which  will  be  chosen  later.  Then 

K  -  C  *  I 

and  Eq.  (All)  becomes 

*  5*  -t'K 


or 


»-f 


CK  +  2(pjg  CM^  co) 


In  order  to  use  the  Poisson  integral  formula,  K"  is  chosen  such  that 

zrr 


I. 


a.  • ; 


cL  u)  =  0 


which  leads  to  K  =  “  <^,g  ,  so  that 


(H  » 1 


=  4>jg(l~2c*i  co) 


(A13) 
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Finally,  by  Poisson's  integral  formula  (see  Reference  15,  p.  254^ 

$  j  =  -  — —  I  (  1  ~  2  Cm  ^  lo')  Jnx.  [l  -  2  ej)  +  ]  ci  w' 

>'<9 

The  function  is  now  determined  up  to  a  fimction  of  ^ say,  and 
possible  eigensolutions  that  are  required  to  match  with  the  outer  solution. 


A. 1.2  Matching  and  Final  Solution  Forms 


The  outer  solution,  Eq.  (A4) ,  when  written  in  terms  of  the  inner 
variables,  becomes 


‘p^'°  -  t/  *  + 


2n  J 


d  _ 


Ve  j.. 


When  this  equation  is  expanded  for  small  S  ,  using  Eq.  (AS),  the  result  is 
0^'®  ~  a„x  +  £U„  i-6<t>,,ci)y  +  6^<p,^(X}V^ 


P  -  /  PCz^)  d  X 


371 


-  £?cf; 


-  Q(.x)  Jivi  /t*  +  J  z{x-x^)d.x^ 

-flO 

~  I  ^  oi6\<r^i£)  s)^ 


15.  Churchill,  R.V.  Complex  Variables  and  Applications,  Second  Edition, 
McGraw-Hill,  New  York,  1960. 


(A14) 


CA15) 


(A16) 
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The  inner  solutionj  Eq.  (A6] ,  when  written  in  terms  of  the  outer 
variables  and  using  the  results  in  Eqs.  (AlO]  and  (A14),  becomes 


~  ^  /  (t  -  2  u>' )  [,1  ~Zi^/S) 


+  (.r/ S')*!  dui^  +  S* 

When  this  equation  is  expanded  for  small  S  ,  the  result  is 


CA17) 


^  ~  Ll„  z:  + 


^  +  [<>«  ^  +  <5 


(r/S)  -f-S^c^gCziJ  ^  o(d*)j-  (AIS) 

Matching  Eqs.  CA16)  and  (A18)  results  in 
&,(£)  »  a 
PCz)  =  an 

to  complete  the  solution  for  the  line  doublet  distribution,  and 
O-^(S)  =  S* 

i2C2)  =  7T  (^x) 


~  <i>i2  cz>  £ I5  +  1  f  (;6,g(£,)  a LtC  cL 
"  2  /  ^  2(x^~x  )d.x^ 


to  complete  the  solution  for  the  line  source  distribution.  Inspection  of 
the  first  term  of  Ggixy  indicates  that  an  additional  term  of  order 
should  have  been  assumed  in  the  inner  expansion  of  Eq.  (A63 .  However,  considera¬ 
tion  of  the  boundary  conditions  quickly  shows  that  this  additional  term  is  a 
function  only  of  X.  which  is  precisely  this  portion  of  G-j  CX)  .  Furthermore, 
no  eigensolutions  are  necessary  to  perform  the  matching. 
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A. 1.3  Pressure  Coefficient 


The  pressure  coefficient  may  be  calculated  from  the  inner  solution 
and  to  second  order  in  the  small  parameters,  i.e.,  retaining  €.  ,  6  and  £<S, 

-  2€6(4>'„Y  *  ^2-  )  - 


^  ^  (<P„  +  CA19) 

On  the  pipe  at  W  =  0  and  TT, 

Cp^  ^  -  2  €  4>\^  +  <1.  €  6  <f>'„  ~  ^  ^  (1  ‘  )  (A20) 

Cp^  =  -Z£  (A213 

In  terms  of  the  physical  quantities,  these  equations  become 

+<<y  v/(i;  ~  j3^u^  {£)  (A22) 

=  -?U.^(.i)-4S  vjli)  -  CA23) 

which  are  the  same  as  had  been  derived  in  the  main  text. 


A. 2  ALTERNATIVE  TWO-DIMENSIONAL  STATIC-PIPE  APPLICATIONS 

The  Static-pipe  response  analyzed  in  Section  2.3  and  Appendix  A.l 
can  be  exploited  to  determine  directly  without  the  necessity  for  a 

separate,  independent  probe  measurement  to  determine  the  constant  of  integra¬ 
tion  in  Eq.  (29)  of  the  main  text.  This  alternative  procedure  would  involve 
a  different  pipe  configuration  from  the  ones  that  were  fabricated.  In  partic¬ 
ular,  a  third  set  of  pressure  measurements  must  be  made  with  each  pipe, 
either  by  a  third  row  of  orifices  or  by  rotation  of  the  pipe. 

This  procedure  can  be  described  in  terms  of  Eqs.  (22)  to  (24)  of 
the  main  text,  which  give  the  pressure  coefficients  on  the  pipe  at  =  0,  Tt/2 
and  TT  ,  called  Cp  ,  Cp  and  Cp  ,  respectively.  These  equations  are 
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u^(£)  ■/■  4  S  zrj  (z;  (A24) 

^f'o  =  f)  =“  -  c£:j>  - /S®u!/c^;  -  4  v/ c:£j  CA25) 


Cp^  5  CpC%,  7tl 


- 2  a,  (£ >  - y3 ® y  -  ^Sv„'cz> 


(A26) 


It  should  be  noted  that  alternatively,  C/»  =  Cp  (i,  ilt/z),  since  the  re¬ 
sults  for  60  *  7t/z  and  57E/2  must  be  equal  by  symmetry.  If  Eqs.  (A243  and 
(A26)  are  added  and  subtracted,  Eqs.  (25),  (26)  and  (28)  of  the  main  text 
are  obtained,  namely 


Cff,  +  Crp^ 

4 


/5  = 


8  ^ 
LX) 


Z 


(A27) 


which  when  solved  for  Lu^  , 
U,  CaTJ  -  - 


(A28) 


and 

^  f  ^ 

(«)  = 

as 


(A29) 


If  Eqs,  (A25)  and  (A27)  are  combined,  (£p 


is  found  directly,  namely 


2  - 


(A30] 


Therefore,  by  making  the  three  sets  of  measurements  of 
and  Cp^  ,  can  be  found  directly  from  Eq.  (A28) ,  "ItJ  LZ)  from  Eq.  (A29) 

and  from  Eq,  (A30).  Finally,  {X  )  can  be  found  from 


V,LX)  =  V^7i7  JU^\^^dvfLZ}/di]/ 


CA31) 
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-/ 

o 

1 


X  <  o 

X  ^  o 
X  >  o 


(A32) 


In  principle,  it  would  only  be  necessary  to  measure  at  a  suffi¬ 
cient  number  of  points  so  that  ti)  and  its  derivative  could  be  determined 
accurately  for  evaluation  of  by  Ec{,  tA31)  at  a  single  point,  say  X^  • 

Once  Vg  iXg^  is  known,  Eq.  (A29)  can  be  integrated  as  in  Eq.  C29)  of  the  main 
text  to  obtain  iX)  over  the  entire  interval  for  which  ■Vj'c£)has  been 
measured.  If  additional  Cp  measurements  are  made,  the  accuracy  of  the  entire 
scheme  can  be  assessed  by  comparing  the  results  for  'tr^  (£)  as  determined  by 
Eq.  CA31)  with  those  found  by  integrating  Eq.  {A29) . 

A. 3  STATIC-PIPE  APPLICATIONS  IN  THREE-DIMENSIONAL  FLOWS 

The  ^'and  ^  components  of  the  three-dimensional  flow  that  is  in¬ 
cident  to  the  pipe  are  shown  in  the  sketch,  where  both  '^q  and  <X) , 
and  so  and  (X'i  are 

arbitrary  functions  of  X  .  The 
relationships  among  these 
quantities  are 

-  ^(xy  uj^cx')  CA33) 

iXy  »  %  iX)  GC;  (A34] 

c  cj^  (X)  {A353 

O  w  9 
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The  analysis  given  in  the  main  text  can  be  adapted  readily  to  this 
configuration.  In  particular,  the  potential  for  the  crossflow  in  Eq.  (5)  is 
replaced  by 


(A36) 

whereupon  Eqs.  (6)  to  (8)  are  replaced  by 

v^cx.r.uj) 

*  C-  R/r)  d  ^%(X)  [iAJ  /d.TC’ 

(A37) 

=  X cx)  R ^ £uj  ~oj^  cx)]  / 

(A38) 

v^^z,r,uj) 

=  a;>  //■  ® 

(A39) 

If  the  analysis  is  followed  through  exactly  as  in  the  main  text,  Eq.  (21)  is 
replaced  by 

<=  -  2o~„(k)  +  ^6  d  £uj^  /dx 

-  To  (.'X') Z,w-oJ,CX.}  2  CA40) 

or,  equivalently,  using  Eqs.  CA33)  and  CA34), 


CpiX,  uf)  =  "  z  u.^  ix)  -  /3^  iX.)  +  4Sd  [.Vg  (.X)  /di 


-  <  C*J  -  C^To  CX  )  C-W  t*J  J 


CA41) 


Measurements  must  be  made  every  90“  around  the  pipe  for  three-dimen¬ 
sional  flows,  either  by  four  lows  of  static  pressure  orifices  or  by  a  lesser 
number  of  rows  and  pipe  rotation.  For  measurements  atw  =  0,  7T  ,  rf/s  and 
3  7Z/2,  it  follows  that 

Cp^  s  Cp(x^o)  «  -  2u,Cx)  +45  v'o(.^)  “4  1  (A42) 

-  G.p{.X.Tt)  -  -  2  U^(X  )  - llI  (.X  )  -  4  <5  vj  (i)  -4  ex'}  (A43) 


64 


AEDC  TR-80  34 


s  CpC%,  *  -  2  a^cx)  - CZ)  +  <5  -  4  ir/ )  (a44) 

-^)  =  -  -/3  -  4-^6tr  C*j  -  (.X){MS'^ 

Adding  all  four  of  these,  and  using  Eq.  CA^S^, 


*  _  /5*ct/c£j  Cf,  *  i-Cfl  +  Cf,^ 

U;  «)  +  - j— - g -  -  <*' 


which  can  be  solved  for  t*^C;£^to  give 
_  ^  \  V  li^ 


A  2 


*  -  ^.  {r-  [r-  P 

Subtracting  from  and  then  Cp^  from  give,  respectively, 


v.  CZ) 


^P,~^P-L 

8S 


CA46) 


CA47) 


In  addition,  the  U^fx),  2^'cxj  and  ^ c  x i  terms  can  be  eliminated  from 
Eqs.  (A42)  to  (A45)  to  give 


(A48) 

(A49) 


*  VA  VX- 


^  8 


CA50) 


With  the  relationships  of  Eqs.  (A47)  to  (A50) ,  then,  if  measurements 
of  all  four  pressures  are  made  along  the  pipe,  and  at  least  one  independent 
measurement  of  either  (2)  or  cD-g,C2)  is  made,  say  by  a  probe,  all  of  the 
/  desired  velocity  components  can  be  determined.  For  example,  if  Xr^  (.2^)  is 

measured  independently,  Eq.  (A48),  which  is  identical  with  the  two-dimensional 
equation,  can  be  integrated  according  to  Eq.  (29)  of  the  main  text  to  give 
(.X) .  Then  Eq.  (ASO)  can  be  solved  for  C.Z  )  ,  which  in  turn  can 

be  solved  for  uT^  (^*,  using  Eq.  CA49),  by 

Jy.  /oix'j  y (A51) 
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''8 

where  sgn  ^  is  defined  in  Eq,  CA32).  Finally,  can  be  found  from 

Eq.  CA35)  whereupon  (.X)  can  be  found  from  Eq.  CA475 .  As  an  alternative, 
of  course,  5  could  be  measured  independently  and  a  similar  procedure 


foUowed  to  determine  the  distributions  of  IX) ,  (Jc)  , 

^  )  and  (X^  CtC  ) . 
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NOMENCLATURE 

pressure  coefficient,  -  fVo  )/<?«> 
airfoil  chord 

functions  in  definitions  of  and  ,  Eqs.  (A9)  and  (A12), 
respectively 

height  of  static  pipe  above  tunnel  centerline.  Figure  2 
constant  in.  definition  of  ,  Eq.  (A12^ 
iterative  relaxation  factor,  Eq.  (30) 

Mach  number 

doublet  strength  in  outer  solution  for  Eq.  (A4) 
static  pressure 

source  strength  in  outer  solution  for  <p  ,  Eq.  (A4) 

free-stream  dynamic  pressure, 

nondimensional  radius,  r/ R 

static  pipe  radius.  Figures  2  and  3 

Reynolds  number  based  on  airfoil  chord 

radial  coordinate.  Figure  3 

function  describing  surface  of  static  pipe 

signum  function,  Eq.  (A32) 

free-stream  velocity 
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ti-jir,  disturbance  velocity  components  in  « ,  4  j  ^  directions, 

respectively 

^Critical  ct  for  which  U^  +  ll  corresponds  to  1 

nondimensional  a  found  from  external-flow  calculation 
three-dimensional  velocity  normal  to  static  pipe,  Eq.  (A355 


nz 

} 

Oi 

/3 

<5 


cL  v.  /dX 


disturbance  velocity  components  in  X,  r  ,  u)  directions, 
respectively,  Figure  3 


rectangular  coordinates  with  origin  at  tunnel  centerline  and 
u.  axis  coincident  with  line  between  junctions  between  plenum 
cnambers  6  and  7  and  16  and  17,  Figures  1  and  2 

coordinate  from  airfoil  leading  edge 

reference  value  of  X  for  determination  of  in  Eq,  C29) 

nondimensional  coordinates  and  ,  respectively 

coordinate  from  static  pipe  centerline  equal  to 
Figures  2  and  3  and  Eq.  (2)  ^ 

static  pipe  coordinate,  Figure  3 

airfoil  angle  of  attack 


V7T 


Ml 


nondimensional  static  pipe  radius,  f?/c 


£  parameter  representing  magnitude  of  model/wall- induced  dis¬ 

turbances  due  to  thickness  or  angle  of  attack 

air  density 
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$ 

4> 


^10  •  ■•  ^t2 

6iJ 


7 


gauge  function  in  outer  solution  for  <P  ,  Eq.  (A4] 
airfoil  thickness -to-chord  ratio 
inner  solution  velocity  potential,  Eq.  CA6) 
terms  in  series  expansion  of  Eq.  CA6) 
function  in  definition  of  ,  Eqs.  (A123  and  (A143 
perturbation  velocity  potential 
functions  in  expansion  of  ,  Eq.  (Al) 
functions  in  expansion  of  Eq.  CAS] 

angular  variable  in  cylindrical  coordinate  system,  Figure  3 

angular  orientation  of  three-dimensional  disturbance  flow 
incident  to  static  pipe.  Appendix  A. 3 

gradient  operator 


SUPERSCRIPTS 

m  model/wall  induced  quantity 

^  static  pipe  induced  quantity 

T  total  quantity 

T,I  total  inner  solution  quantity 

J.O  total  oufer  solution  quantity 

ft), (3^  iterative  step  number 

calculated  distribution  before  relaxation  to  take  second 
iterative  step 
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^  nondiraensionalization  of  velocity  through  division  by 

"  nondiraensionalization  of  length  through  division  by  ^ 

SUBSCRIPTS 

m  measured  quantity 

I 

O  quantity  at  static  pipe  centerline 

^2,3, 4  quantity  at  bottom,  top,  left-hand  side  and  right-hand  side, 
respectively,  of  static  pipe;  Ui  =  o,7r  ^jr/s  and  371/2  , 
respectively,  in  Figure  3 

<»  free- stream  quantity 
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